Experimental observation of second-harmonic generation and diffusion inside random 

media 
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We have experimentally measured the distribution of the second-harmonic in- 
tensity that is generated inside a highly-scattering slab of porous gallium phos- 
phide. Two complementary techniques for determining the distribution are 
used. First, the spatial distribution of second-harmonic light intensity at the 
side of a cleaved slab has been recorded. Second, the total second-harmonic 
radiation at each side of the slab has been measured for several samples at vari- 
ous wavelengths. By combining these measurements with a diffusion model for 
second-harmonic generation that incorporates extrapolated boundary condi- 
tions, we present a consistent picture of the distribution of the second-harmonic 
intensity inside the slab. We find that the ratio i2u/Lc of the mean free path at 
the second-harmonic frequency to the coherence length, which was suggested 
by some earlier calculations, cannot describe the second-harmonic yield in our 
samples. For describing the total second-harmonic yield, our experiments show 
that the scattering parameter at the fundamental frequency /cic^^ioj is the most 
relevant parameter in our type of samples. 
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1. Introduction 

Optical second-harmonic generation inside inhomogc- 
neous media has attracted much interest for biologi- 
cal applications such as high-contrast microscopy [1| 
and for photonic applications such as high-efficiency fre- 
quency conversion in granular nonlinear dielectrics 0. 
The second-harmonic signal radiated from a random 
medium is also an informative probe of the fundamen- 
tal and the second-harmonic intensity distribution inside 
the medium. Second-harmonic signals provide valuable 
information about multiple-scattering processes inside a 
random medium, which is not easily accessible by other 
methods. Previously, these signals have been used to 
study the angular, spatial, and temporal correlations in- 
side scattering media 0, |3j Hi ■ 

Currently available first-principles theories of optical 
second-harmonic generation in multiple scattering media 
consider only the case of the mean free path being much 
larger than the wavelength, far from the Anderson local- 
ization regime [1, 0|- Many articles have discussed in- 
terference effects in multiple-scattering nonlinear media 
such as the effect of weak localization or the enhanced 
forward scattering [1, [^, [l^, [U . 

The efficiency of the nonlinear conversion in granu- 
lar dielectrics is related to the minimized role of phase- 
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matching in such a medium. Phase-matching is essential 
for the efficiency of three- wave mixing processes in non- 
linear materials. The optical dispersion in the nonlinear 
material results in a phase-mismatch between the funda- 
mental and the higher-harmonic propag ating hght [ij]. 

There are two conventional methods to overcome the 
phase mismatch. The first method makes use of the 
fact that in birefringent and nonlinear materials, it is 
possible to fulfill the phase-matching condition by ad- 
justing the crystal orientation with respect to the beam. 
The second method, called quasi-phase-matching, uses 
a periodically-polled polarizability of the material. The 
nonlinearity of the crystal is modulated such that the ex- 
tra momentum, which results from the phase-mismatch, 
can be transferred to the crystal. In both of these meth- 
ods, the total conversion yield increases quadratically 
with the length of light path inside the medium. 

Recently, a novel method, called random-quasi-phase- 
matching, has been introduced in which the crystal ori- 
entation is randomly varied along the beam path by us- 
ing a powder 0]. As suggested from theory and con- 
firmed by the experiment, for this new method, the accu- 
mulated second-harmonic energy increases linearly with 
the sample thickness. By using this method, a rela- 
tively high second-harmonic signal can be extracted from 
many semiconductor powders, which were otherwise (in 
their bulk form) useless for conventional methods. Using 
random-quasi-phase-matching, second-harmonic genera- 
tion can be obtained for a larger bandwidth and accep- 
tance angle than conventional methods. 
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A high second-harmonic yield has also been observed 
by Tirinyanu et al. for strongly-scattering porous- 
GaP |l3| . Mernikov and coworkers have shown that the 
second-harmonic signal in the specular-reflection direc- 
tion is enhanced by orders of magnitude after anodically- 
etching the single-crystalline GaP wafer [3, EE] ■ How- 
ever, no comparison with a theoretical model had been 
provided for these observations. The distribution of 
second-harmonic light inside such a medium has also not 
been investigated experimentally. 

In this article, we use two complementary techniques to 
determine the distribution of the second-harmonic inten- 
sity inside a multiple-scattering slab. In the first exper- 
iment, the spatial distribution of second-harmonic light 
intensity at the side of a cleaved slab is measured. In the 
second experiment, the total second-harmonic radiation 
at each side of the slab is measured. On the macroscopic 
level, our experimental results confirm the predictions 
of the diffusion theory for distribution of the second- 
harmonic intensity in a multiple-scattering medium. On 
the microscopic level, however, the only available theo- 
retical model by Kravtsov et al. @, which is based on 
the diffusion theory, is unable to describe our experi- 
mental results. Their model suggests the ratio £2uj/Lc 
of the transport mean free path at the second-harmonic 
frequency to the coherence length as the universal pa- 
rameter for the conversion yield. By performing mea- 
surements at a range of scattering strengths and frequen- 
cies we find that their model does not apply to our type 
of samples. Instead, our measurements show that the 
second-harmonic yield follows a consistent dependence 
on the transport mean free path at the fundamental fre- 
quency ilu,. 

In the following sections we briefly review the available 
theories of second-harmonic generation in random me- 
dia. We present our derivation of the intensity distribu- 
tion from a diffusion model with the extrapolated bound- 
ary conditions. In the experimental section we describe 
two separate measurement techniques which probe the 
distribution of second-harmonic light generation inside a 
highly-scattering sample. Results of our experiments are 
then compared with predictions of the diffusion model. 

2. Theory 

A. Optical nonlinearity in opaque material 

An opaque medium may also be optically nonlinear. This 
nonlinearity can be an intrinsic property of the bulk ma- 
terial or a result of the enormous interfacial area present 
in porous objects. The second-order nonlinearity is ab- 
sent in many non-crystalline materials or crystal struc- 
tures because of the presence of inversion symmetry. Rel- 
atively large non-linearities may arise at the interfaces 
of these materials with others or with vacuum due to 
symmetry-breaking at the interface. The scientiflc un- 
derstanding of optical nonlinear processes in strongly- 
scattering materials is still in the preliminary phase. 
Some models have been developed [g, 0] based on the 



diffusion approximation, in which interference effects are 
assumed to be averaged out and the sample size L is 
much larger than the transport mean free path. 

In these diffusion models the incident wave at the 
fundamental frequency u experiences several scattering 
events before leaving the random medium. Light at 
the second-harmonic frequency is generated during the 
multiple-scattering process. The propagation direction of 
the second-harmonic light is scrambled within one trans- 
port mean free path £2^, thus becoming an isotropic 
source of diffusive photons at the second-harmonic fre- 
quency. The effect of phase-mismatch between the funda- 
mental and the second-harmonic light is negligible when 
the transport mean free path is much smaller than the 
coherence length Lc{uj) = 7r/|2fci^ — A:2w|, where ki^ and 
k2uj are the wave-vector magnitudes at the fundamen- 
tal and the second-harmonic frequencies. Therefore, in a 
random medium that consists of grains showing nonlin- 
ear response, the effect of constructive interference can be 
overcome by selecting grain sizes to be smaller than the 
coherence length. It has been experimentally shown that 
the second-harmonic yield from equal amount of material 
increases with grain size until the grain size approaches 
the coherence length [2] . 

Overcoming the destructive interference due to phase- 
mismatch is also possible by introducing scatterers inside 
a homogeneous nonlinear crystal. In such a medium the 
fundamental and second-harmonic waves scatter differ- 
ently, therefore the destructive interference of the oth- 
erwise copropagating waves does not occur, providing 

B. Review of existing theories 

Second-harmonic generation in random media has been 
theoretically modeled for three different systems in two 
reports. Here we briefly introduce and compare these 
models as they will be needed for comparison with our 
experimental results. The important prediction of these 
models is the magnitude of the mesoscopic jll] second- 
harmonic conversion rate T which is defined as the 
second-harmonic intensity per volume, generated inside 
the scattering medium, per unit of fundamental energy 
density squared . 

In the first report Kravtsov et al. Q have considered 
two different systems. In the first system, scattering is 
introduced in a slab of nonlinear dielectric by adding a 
small volume fraction of point-like scatterers. The mi- 
croscopic nonlinear polarizability is assumed to be unaf- 
fected by the presence of scatterers, thus constant across 
the sample. We refer to this model as "homogeneous 
nonlinear background model" . The medium can be opti- 
cally dispersive. An = n2u) — «iw 7^ 0. When Art < 0, the 
conversion rate T is independent of scattering parameters 
kiujU-iu and k2u,£2ui- When An > 0, F is found to depend 
on ki^£i^ as 

r = Asou;n2^nUx^'^r arccot[-^^%— ], (1) 

^2 + itzi^'^ 
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where A is a dimensionless prefactor, x^^^ is the norm 
of the nonlinear polarizability of the bulk medium av- 
eraged over all directions, Lc is the coherence length, 
and ^lui and are mean free paths at fundamental and 
second-harmonic frequencies. All parameters in Eq. ([1]) 
can be frequency dependent due to the optical dispersion 
of the material. In Eq. ([1]), x^^-* and Lc are microscopic 
properties of the bulk material while ti^ and (,2u] are 
mesoscopic properties that depend on the porous struc- 
ture of the material. Equation ([1]) predicts that as the 
mean free path at second-harmonic frequency becomes 
small relative to the coherence length, the conversion rate 
should increase. The derivation of this result assumes 
weak multiple-scattering, fcic^^it^ ^ 1, A;2cj^2ij S> 1. 

The second system that Kravtsov et al. [6] have con- 
sidered is a dense powder of nonlinear grains in which 
the grain size is larger than the wavelength and much 
smaller than the mean free path. This system is similar 
to samples analyzed in reference [2]. We call this model 
"nonlinear-powder model" . 

In the second report, Makeev and Skipetrov 0] have 
introduced a third system. They have modeled a sus- 
pension of colloidal particles. We refer to this model as 
"nonlinear colloidal suspension model". In their model, 
the conversion centers are the same as the scatterers and 
the background medium is linear. They found that, for 
these suspensions, the second-harmonic intensity divided 
by the number of scatterers shows no explicit dependence 
on the multiple-scattering properties of the suspension. 
In their model, the second-harmonic conversion rate is 
given by the following simple product: 



r = BpV2u,^2Lj, 



(2) 



where p is the concentration of colloidal particles, S2aj 
is the total second-harmonic cross-section of an isolated 
particle, V2uj is the energy velocity at second-harmonic 
frequency, and B is a dimensionless number. 

In all three models, the nonlinearity is assumed to be 
so low that the distribution of fundamental diffusive pho- 
tons is not affected, therefore the energy density at fun- 
damental frequency Uiu,{z,t) can be calculated by solv- 
ing the diffusion equation for an optically-linear multiple- 
scattering slab. 

C. Macroscopic distribution of the second-harmonic in- 
tensity inside an opaque slab 

Following the literature [6l, 't'I , we assume the source dis- 
tribution of the second-harmonic diffusive photons inside 
the opaque material S'2w(r, i) to be equal to the conver- 
sion rate F times the square of the energy density of the 
fundamental light Uiui{r,t): 



^2c.(r,i) =^^7L(r,^)• 



(3) 



The exact value of F depends on the model and the sam- 
ple type [Eq. dJ) sets an example]. Being a mesoscopic 
quantity, F will not affect the macroscopic intensity dis- 
tribution except for changing an overall prefactor. 



For a slab of finite thickness, illuminated with a con- 
tinuous plane wave of intensity Iq, the following set of 
diffusion equations can be written for the distribution of 
fundamental and second-harmonic diffusive-photon den- 
sities, 



D2 



' dz^ 
d^U2M) 
' dz^ 



= -ruUz), 



(4) 
(5) 



where Di^ = vi^ii^/S and D2i^ = i'2w4w/3 are the dif- 
fusion constants at fundamental at second-harmonic fre- 
quencies, vii^ and V2uj are the energy velocities in the 
medium. The source term Si^j^z) denotes the distribu- 
tion of the diffusive light source at the fundamental fre- 
quency. These equations must be solved together with 
the following boundary conditions [13], 
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where ^(^, are the extrapolation lengths at the inci- 
dence interface of the slab, and ^('^, C,2uj are the extrapo- 
lation lengths at the opposite interface of the slab at the 
corresponding frequencies. 

There are various ways of formulating the source distri- 
bution at fundamental frequency Siuj{z). A phenomeno- 
logical and perhaps the most practical method [3, 
is to set the source at one transport mean free path 
from the incident interface, inside the slab: Siaj{z) = 
{viljIo/^1ui)S{z/£iui — 1). A natural extension [l^ 
to the previous description considers a source term 
that exponentially decreases with depth, with a de- 
cay length equal to the mean free path: Si^{z) = 
(vi^Iq/Ii^) exp{—z/£i^). In the following calculation we 
use the latter description. 

By integrating the diffusion equation ^ over an ex- 
ponential source and applying boundary conditions ([6|) 
and ([7]) we derive, for inside of the slab. 



+ CL){L + CL - ^) 



o exp(--^)], 
(10) 

where the thickness L of the slab is considered much 
larger than the mean free path, and thus terms such as 
exp(— L/^oj) are neglected. 

The second-harmonic energy density is given by in- 
serting the result of Eq. (fTU]) in diffusion equation ([5]) 
followed by a double integration. The total generated 
second-harmonic radiation that is propagating in the 2t: 
spherical angle around the incidence direction is referred 
to as the forward radiation and denoted by T2uj ■ The total 
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radiation that is propagating in the opposite (backward) 
direction, is denoted by i?2i^- For an slab illuminated by 
a plane wave, these two quantities are given by 



^2u — H -1 U=0 : 

V2ui dz 

_ D2u> dU2oj{z) . 

^2uj — ; Z=i 

V2ui dz 



(11) 

(12) 



Inserting the fundamental photon-density distribu- 
tion pop in diffusion equation ([5]) , and integrating yields 
the second-harmonic photon-density distribution inside 
the slab, 



U2M) = {C,L + z)R2 



Ufjz2)dz2dzi. (13) 



Boundary condition ([S]) has been applied. Applying the 
boundary condition ^ yields R2uj- Expressions for U2uj, 
R2uj, and T2U1 as functions of T, /q, L, -^ic^, Ciw, and (2^^ 
were obtained by combining Eqs. (fTO|) . (fTTj) . and (fT3l) 
with boundary condition ^ in Mathematica. The closed 
forms of the answers contain many terms and are not nec- 
essary to present here for the purpose of the discussions 
in this paper. In completion of earlier calculations 7], 
we have included the extrapolated boundary conditions 
and exponentially decaying source term for fundamental 
light. 

We have compared the internal distribution of second- 
harmonic intensity in the limiting case of vanishing ex- 
trapolation length — C2UJ — with the corresponding 
result stated in Ref. 7] and found that their approach 
works well for the case of optically-thick slabs with zero 
extrapolation length, i.e. absorbing boundaries. The ef- 
fect of finite extrapolation length is not negligible, re- 
gardless of the slab thickness. 

Here we present the answer up to the first order in 
£i^/L. The first two non-zero orders sufhcc for most of 
our discussions. 



R2 
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As has been mentioned earlier in this section, the meso- 
scopic conversion rate F shows up only in prefactors in 
Eqs. and (|15p and does not affect the macroscopic 
distribution of radiation inside and around the slab. Our 
derivation shows that the total second-harmonic intensity 
is enhanced with increasing extrapolation ratio Ciujl^i^- 
Physically, this enhancement is caused by the increased 
trapping of fundamental light inside the diffusing medium 
due to the interfaces. 



Due to the variability of several parameters from sam- 
ple to sample, it turned out to be useful to experimentally 
determine the ratio 



R2U1 
T2L0 



(16) 



In theory, the dimension-less quantity "q is independent 
of the extrapolation ratio. Up to the the first order in 
(^ii^/L), 77 is also independent of the conversion rate 
F and the incident intensity /q. As reported in earlier 
works ji], for optically-thick slabs, ratio 77 converges to a 
constant, i.e. lim 77 — > 3 when L/ —f 00. 

The most appropriate (currently available) theoretical 
model for our type of samples is homogeneous nonlin- 
ear background model described in Sec. 12 B[ in which 
the conversion rate is given by Eq. ([T]). The intensity 
of the second-harmonic radiation is dependent on sev- 
eral parameters. To compare any experimental result 
with any multiple-scattering theory, it is advantageous to 
separate the dependencies induced by multiple-scattering 
(mesoscopic) from variations that are caused by intrin- 
sic material properties (microscopic) and from geomet- 
rical specifications of the sample (macroscopic). There- 
fore, for comparing the experimental results with theory, 
we divide the measured second-harmonic backward radi- 
ation i?2a) by the sample thickness, the incident inten- 
sity squared and all prefactors of the arctan function in 
Eq. Il]) , which are material properties of GaP. We define 
the normalized yield as: 
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R2 



ujn2u>n^ 



?Lir 



(17) 



which is useful for comparing our data with the theories. 

Note that in all of the above calculations, the absorp- 
tion of fundamental and second-harmonic light is ne- 
glected. This assumption is justified when the absorption 
length of the medium is much larger than LP' /iiui, which 
a justified assumption for our type of samples. 

For a pulsed light-source, generally, the stationary cal- 
culation is valid if tq 3> t^i , where tq is the pulse duration 
and td is the Thouless time, defined as = LF' /Di^. 
Using analytical calculations, which will be presented 
elsewhere, we have calculated the effect of pulse dura- 
tion on second-harmonic intensity distribution. Although 
the non-stationary calculations are analytic, they contain 
long and complicated mathematical expressions that do 
not provide any intuition about the physical outcome. 
We found that for our specific experimental conditions, 
the difference between non-stationary and stationary so- 
lutions is much less than our experimental error, there- 
fore we have presented our results only in comparison 
with the stationary solutions. 

3. Samples and setup 

A. Samples 

The porous-GaP samples were etched from commer- 
cially available single crystal n type gallium phosphide 
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Table 1. Summary of specifications of the samples that 
are analyzed in this report. T4 is the applied electromotive 
force during the etching process. The measured value of 
mean free paths £21^ and £1^1 are presented at Aq = 0.65 
and 1.3 /im, correspondingly. L is the thickness of the 
porous region. 

1 35.5 ±1 0.27 ± 0.04 2.49 ±0.29 14.7 

2 43.7 ±1 0.26 ± 0.03 1.85 ±0.29 14.7 

3 23.5 ±1 0.26 ±0.04 2.18 ±0.24 14.7 

4 32 ±1 0.83 ±0.06 2.50 ±0.28 22.5 

5 26 ±1 0.60 ±0.05 2.22 ±0.18 11.2 

6 83 ±3 1.10 ±0.08 4.71 ±0.37 10.0 

7 126 ±4 1.24 ±0.08 3.38 ±0.48 15.0 

8 116 ±3 1.58 ±0.11 5.75 ±0.60 10.0 



wafers [20l|. These wafers were doped with donor atoms 
with a density range of 10^^ cm~^. The fabrication pro- 
cedure is as follows. One side of a 9 x 9 mm^ piece of GaP 
wafer with a thickness of 500 /im is anodically etched in 
sulfuric acid under a controlled electromotive force. The 
electric current is recorded in time. The total charge 
that is transported in the process is proportional to the 
amount of removed material, allowing for the calcula- 
tion of the porosity for a known thickness of the etched 
layer. Different post-etching processes can be used to 
remove the top layer, in order to improve optical char- 
acterization, or to increase the pore-size. We used eight 
samples with different values of mean free path. These 
samples were prev iously fabricated for other optical ex- 
periments [mill,!!!. 

A summary of the scattering characterization proper- 
ties of the investigated samples is shown in table 1 . The 
thicknesses were determined from scanning electron mi- 
croscope images. The mean free path values were deter- 
mined by the standard technique of measuring the linear 
total transmission. As a representative, the mean free 
path of sample 2 has been determined over the spectral 
range of 0.6 to 1.6 /im (Fig. [T]) by coupling a spectrom- 
eter to an integrating sphere setup. Dependence of the 
mean free path on wavelength is largely due to the nano- 
structure of the samples. It also varies from sample to 
sample. The magnitude of the nonlinear susceptibility 
of porous GaP for some lattice orientations is compara- 
ble to that of a commercial second-harmonic crystal such 
as potassium dihydrogen phosphate (KDP). However, as 
GaP is not birefringent, it has no application for conven- 
tional methods of second-harmonic generation. 

B. Setup for measuring total radiation 

Fig. [3 shows a diagram of the integrating-sphere setup 
that was used for measuring the total second-harmonic 
light that is radiated in the forward or backward direc- 
tion from a porous-GaP slab. Several gold mirrors guide 
the beam vertically into the entrance of a Ti02 coated in- 
tegrating sphere. The sample can be laid down at the en- 
trance for transmission measurements or placed beneath 
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Fig. 1. The mean free path of a Porous-GaP sample is 
wavelength-dependent. As a representative, we plot the 
mean free path as a function of vacuum wavelength for 
sample 2. Similar curves are also obtained for all samples 
using total-transmission measurements. 



the sphere for reflection measurements. Two photodiodes 
are positioned at two exit holes of the integrating sphere. 
At one exit-hole, an amplified silicon detector (PDA55a - 
Thorlabs) connected to an oscilloscope (DL9040L - Yoko- 
gawa) measures the second-harmonic signal after the fun- 
damental infrared light is filtered by a cold glass filter 
(KG5). At the other exit-hole, a germanium detector 
(DETIOA - Thorlabs) connected to another channel of 
the same oscilloscope, measures the scattered fundamen- 
tal light. In this way, both fundamental and the second- 
harmonic signals can be measured simultaneously. 

For a second-harmonic signal to be detectable by nor- 
mal silicon detectors, an incident intensity of several 
GW/cm^ is essential. A long pulse at this intensity will 
damage GaP, therefore pulses shorter than nanoseconds 
with low repetition rates are needed. The light source 
we used was a traveling-wave collinear optical paramet- 
ric amplifier of super-fluorescence (TOPAS - Light Con- 
version) pumped with a sub-picosecond laser (Hurricane 
- Spectra Physics) at 800 nm and a repetition rate of 1 
kHz. The output central wavelength could be tuned con- 
tinuously from 0.55 to 2.2 /im with a bandwidth of 5 nm. 
The pulse was transform-limited with a duration of 150 
fs. The undesired frequencies are filtered out from the 
output beam by using a polarizer and a low-pass filter 
(RG850 - Thorlabs). 

For each sample, the radiated second harmonic light 
in the forward and backward directions was measured in 
a fundamental- wavelength range of 1.2 to 1.6 /im. Be- 
fore each measurement, the incident power was checked 
by a pyroelectric head and a power meter (Ophir). The 
transmittance of filters and responsivities and the linear- 
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Fig. 2. Experimental setup for measuring the total 
second-harmonic radiation in forward and backward di- 
rections. The porous-GaP slab is illuminated with a par- 
allel beam of 150 femtosecond infrared pulses. The inte- 
grating sphere collects the fundamental and the second- 
harmonic light radiated in all directions from one side of 
the sample. The sample position is either on the top of 
the sphere for transmission measurements or attached be- 
neath the sphere for reflection measurements. A Silicon 
photodetector (PD-Si) behind a cold glass filter measures 
the second-harmonic signal in the visible range. In par- 
allel, a Germanium photodetector (PD-Ge) measures the 
transmitted fundamental light at infrared range. 



ity of detectors were carefully taken into account before 
extracting the second-harmonic yield. 

Special attention is needed to correct for the transmit- 
tance of the non-porous substrate. The total transmit- 
tance of the porous-bulk interface depends on the refrac- 
tive index mismatch at the interfaces and the directional 
distribution of outgoing diffuse light known as the "es- 
cape function." We used the escape function presented 
in [23|j with an effective refractive index of 1.6 ± 0.2 for 
the porous region. This effective refractive index was 
measured before based on the filling fractions and es- 
cape function measurements from reference Using 
the Fresnel equations for transmittance at substrate-air 
interface, the total angularly-averaged transmittance of 
the GaP substrate for these samples has been calculated 
to be 0.28 ± 0.05, for which the magnitude of the error is 
mainly due to the uncertainty in the effective refractive 
index of the porous medium. 

C. Setup for measuring effusion function 

In order to investigate the second-harmonic intensity dis- 
tribution inside a porous GaP sponge in greater detail, 
the second-harmonic radiation from the side of a cleaved 
porous GaP slab was imaged under a microscope while 
the slab was illuminated by a parallel beam of infrared 
pulses at normal incidence. We refer to the diffusive en- 
ergy flux at the interfaces of the random medium, which 
can be measured with our method, as the "effusion func- 
tions^ . The measurement of the effusion function is sensi- 
tive to details of the distribution of the second-harmonic 
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Fig. 3. Experimental setup for effusion microscopy mea- 
surements. The sample consists of a thin porous and 
highly-scattering layer laying on top of a transparent sub- 
strate. The CCD-camera images the second-harmonic 
signal that is radiated from the narrow cross section of 
the porous part of the sample while it is illuminated by 
a parallel beam of infrared pulses. 

source inside the slab. In order to test this sensitivity, 
we calculated the second-harmonic effusion function for 
a half-slab numerically and compared it with the bulk 
second-harmonic energy distribution in a slab, predicted 
by Eq. We observed that the second- harmonic ef- 

fusion pattern approximates the bulk second-harmonic 
intensity distribution well. This result suggests that one 
can get an immediate qualitative sense of the second- 
harmonic energy distribution from the observation of the 
second-harmonic effusion function. 

The effusion microscopy setup is shown in Fig. [31 The 
incident beam in this setup is a Gaussian infrared beam 
of 3 mm in diameter with its center aimed toward the 
edge of the sample. A cold glass filter in front of the 
CCD camera blocks the fundamental light so that only 
the visible second-harmonic light is captured. 

4. Results 

A. Dependence on incident power 

At an incident wavelength of 1.2 ^m, where the efficiency 
is the highest, the dependence of second-harmonic yield 
on the incident power is measured for several samples in 
the forward direction. The fundamental power is linearly 
proportional to voltage Vq of the silicon detector. The 
second-harmonic power is linearly proportional to voltage 
V2 of the germanium detector. The measured relation 
between output voltages is plotted in Fig. [D From this 
plot a consistent power-law dependence between incident 
and second-harmonic powers is evident for these samples. 
The result of fitting shows a power law P2uj oc Pq with 
a — 1.87 ± 0.03. In a second-order nonlinear process the 
second-harmonic intensity is proportional to the square 
of the incident intensity. The observed deviation of the 
experimentally measured a from 2.0 may be a sign of 
nonlinear (three photon) absorption inside the porous- 
GaP. 
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Fig. 4. The backward-radiated second-harmonic signal 
V2, detected by the sihcon detector, is scaled by the sam- 
ple thickness and plotted for three of the samples ver- 
sus the signal detected by Germanium detector Vq in bi- 
logarithmic scales. The second-harmonic power is pro- 
portional to V2 and shows power-law dependence (ex- 
ponent = 1.87 ± 0.03) on the incident power, which is 
proportional to Vo- Numbers in the legend indicate to 
sample tags as introduced in table 1. 



B. Effusion function at second-harmonic frequency 



The microscopy setup described in Sec. 13 Cl has been used 
to capture the effusion function at the second-harmonic 
frequency. As the incident beam size is much larger than 
the slab thickness, the effusion pattern is laterally (paral- 
lel to the substrate) invariant. Roughness of the section 
due to its porosity and cleaving-defects, can cause fluctu- 
ations in the observed intensity. The intensity distribu- 
tion is averaged in the lateral direction. The experimen- 
tal result is compared with our numeric calculation of the 
diffusion model in Fig. O The bulk distribution of the 
second-harmonic intensity, is also plotted for comparison. 

The measured second-harmonic intensity distribution 
qualitatively agrees with the diffusion model. The inten- 
sity has a maximum around one third of the slab thick- 
ness, as was predicted theoretically in Q. The largest 
deviation is observed near edges, which we attribute to 
the stray second-harmonic light which is leaving the other 
interfaces of the slab behind the imaging plane. 

C. Second-harmonic intensity distribution in far-field 



The measured value of 77, defined by Eq. has been 

plotted versus the optical thickness of samples at vari- 
ous fundamental wavelengths in Fig. [51 The calculated 77 
in the framework of stationary diffusion model, accord- 
ing to Eq. is also plotted for comparison. We find 
a good agreement between measurement and theory for 
most of samples. Large error bars are mainly due to the 
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Fig. 5. Inset: Micrograph of the second-harmonic effu- 
sion intensity from a section of the porous-GaP slab while 
it is illuminated with a parallel beam of infrared laser- 
pulses from left. Brighter regions indicate higher effusion 
of second-harmonic light. Outset: The measured second- 
harmonic intensity is averaged parallel to the substrate 
and its peak is normalized to one. The result of the 
experiment (symbols) is plotted versus the position in- 
side the sample and is compared with the prediction of 
the stationary diffusion model (solid line), found from 
the numeric calculation with no adjustable fitting pa- 
rameters. Our theoretical value for distribution of the 
second-harmonic intensity in the bulk, which is presented 
in Sec. 2C, is shown by the dotted line. The theoretical 
fundamental-frequency intensity distribution inside the 
slab (dashed line) is plotted for comparison. 



uncertainty in the effective refractive index of the porous 
medium. 

D. Second-harmonic yield 

At several incident frequencies, we have measured the 
total second-harmonic intensity radiated from the sam- 
ple. Extracting the absolute value of the conversion rate 
from the second-harmonic intensity measurements is lim- 
ited by a large systematic shift which arises from uncer- 
tainties in the response function of the setup. However, 
this systematic shift does not affect the relative values. 
Therefore for testing theoretical models, we consider rel- 
ative trends rather than the absolute values. 

To check whether the homogeneous nonlinear back- 
ground model of Kravtsov et al. can describe our 
data, the normalized yield 7 defined by Eq. ([T7]) has been 
plotted versus the ratio iiuji^c- This plot is shown in 
Fig. [7{a) for all samples and measured wavelengths. For 
this plot the optical dispersion relation and the nonlin- 
ear polarizability of bulk GaP has been taken from mea- 
surements of ^] and calculations of [2^, respectively. 
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Fig. 6. The ratio rj between total second- harmonic hght 
measured in the backward and the forward direction is 
plotted versus the optical thickness L/ii^^ for various 
wavelengths and samples. The stationary diffusion pre- 
diction from Sec. 2C is plotted as a solid hue. We find 
a good agreement between theory and measurements. 
Numbers 1-8 in the legend correspond to the sample num- 
bers introduced in table 1. 



Our measurements show a trend opposite to the homo- 
geneous background model [H]. For an individual sample, 
the yield increases with increasing mean free path relative 
to the coherence length. The incremental trend is differ- 
ent from sample to sample which indicates that £2uj/ Lc is 
not the universal parameter for describing the conversion 
rate in our kind of samples. 

For the normalized yield plotted versus the scattering 
parameter at the fundamental frequency fcii^^i^, most of 
the measurement points for all the samples are close to a 
single curve. The consistent trend of increasing yield with 
decreasing kii^li^ occurs both when comparing different 
wavelengths in a single sample and when comparing dif- 
ferent samples at the same wavelength. We have fitted 
the data to a power law function and found an average 
exponent of —2.0 ± 0.3. 

5. Conclusion 

Using a microscopy technique, we have measured the 
second-harmonic intensity at the side of a cleaved slab 
during its illumination with a Gaussian beam. We 
have observed that the internal distribution of second- 
harmonic intensity predicted by the diffusion model 
agrees with the experiment. The total intensity radiated 
in the backward direction and the backward-forward ra- 
tio of the second-harmonic intensity has also been mea- 
sured for a number of samples. The measured backward- 
forward ratios show good agreement with the results of 
the diffusion model. 
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Fig. 7. (a) The second-harmonic normalized yield as de- 
fined in Eq. (17) is plotted versus the ratio of the mean 
free path at second harmonic frequency to the coherence 
length i2L,j/ Lc for various wavelengths and samples. The 
normalized yield is defined as the total second-harmonic 
intensity generated in the backward direction divided by 
the square of incident intensity and the thickness of the 
slab and normalized for frequency dependent material 
properties of GaP. The dotted line shows the value cal- 
culated from the theoretical model of Kravtsov et al. @ 
plotted for comparison. No agreement has been found 
between their theory and our measurements. Numbers 
1-8 in the legend correspond to the sample numbers in- 
troduced in table 1. (b) The same data of (a) is plotted 
versus the scattering strength at the fundamental fre- 
quency. The overall trend can be described by a power 
law relation, 7 oc {ki^ii^Y , (3 = —2.0 ± 0.3 which is 
shown by the dashed line. 
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For describing our distribution measurements, we have 
presented a diffusion model for second-harmonic inten- 
sity distribution in strongly-scattering nonlinear mate- 
rial in which the diffused fundamental intensity is con- 
verted into its second-harmonic via the process of de- 
generate two-photon mixing. The second-harmonic light 
also diffuse in the scattering medium. The internal dis- 
tribution and the outgoing intensity in forward and back- 
ward directions of a slab are derived based on the diffu- 
sion equation for light and extrapolated boundary con- 
ditions. As has been reported previously [1,01, diffusion 
theory predicts that in the slab geometry, the total gen- 
erated second-harmonic intensity increases linearly with 
the thickness of the slab and depends quadratically on 
the extrapolation ratio at the fundamental frequency. 
As a result of the diffusion model, the ratio 77 between 
second-harmonic radiated intensity in backward and for- 
ward directions is independent of the mesoscopic conver- 
sion efficiency and is only given by the slab-geometry and 
boundary conditions. Therefore, rj provides a useful test 
for the diffusion model, irrespective of the microscopic 
details. For a sufficiently thick slab, r/ approaches 3 as 
predicted by the diffusion theory. This result has been 
confirmed by our measurements. 

Although current diffusion models describe well the 
distribution of the second-harmonic intensity, the overall 
conversion rate cannot be described by the available the- 
ories. We have found a consistent power law dependence 
between the second-harmonic yield and the scattering pa- 
rameter at the fundamental frequency, 7 oc (ki^uiiuj)^ , 
with an exponent of /3 = —2.0 ± 0.3. 
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